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We show that the magnetic anisotropy in spinel-structure CoC^CU thin films exhibits a strain 
dependence in which compressive strain induces an out-of-plane magnetic easy axis and tensile strain 
an in-plane easy axis, exactly opposite to the behavior reported for the related compound CoFe204- 
We use density functional theory calculations within the LSDA+U approximation to reproduce and 
explain the observed behavior. Using second-order perturbation theory, we analyse the anisotropy 
tensor of the Co 2+ ions in both octahedral and tetrahedral coordination, allowing us to extend our 
results to spinels with general arrangements of Co 2+ ions. 

PACS numbers: 75.70.Ak, 75.80.+q, 75.30.Gw, 71.15.Mb 


Thin films with out-of-plane spontaneous magne¬ 
tization, showing so-called perpendicular magnetic 
anisotropy, are of great interest for applications such 
as high density magnetic memories with fast switch¬ 
ing 1 . Perpendicular magnetic anisotropy is also needed 
for spintronic applications, for example in magnetic tun¬ 
nel junctions, and for low energy current-driven domain 
wall motionSHS. 

In a system of finite size, magnetic anisotropy is deter¬ 
mined by the balance of magnetocrystalline and shape 
anisotropy. The former is a bulk property originating 
from spin-orbit interaction, while the latter originates 
from magnetic dipole-dipole interaction and depends on 
the geometry of the sample. For the case of thin films, 
where shape anisotropy always favors in-plane magne¬ 
tization, it is interesting to understand how the strain 
affects the magnetocrystalline part. Indeed, the way in 
which epitaxial strain can affect the magnetocrystalline 
anisotropy has been the subject of intensive investiga¬ 
tions 6 13 . A prominent example of such investigations is 
the spinel compound C oFe 2 0 4 , which has a strong mag¬ 
netostriction coefficient 7 ® 11 In CoFe 2 0 4 thin films, 
it was shown that changing the sign of the strain leads 
from cooperation to competition of shape and magne¬ 
tocrystalline anisotropies 7 . 

Compounds with spinel (MgAl20 4 -type) structure 
have chemical formula AB^X^ where A and B are cations 
and X represents the ligand anion (usually O, S or Se). 
In this structure, the cations occupy either the tetrahe- 
drally coordinated 15 (T) site or the octahedrally coor¬ 
dinated (O) site. There are twice as many O sites as 
T sites. In a normal spinel, A and B cations occupy T 
and O sites, respectively, while in an inverse spinel, half 
of the B cations occupy the T sites and the remaining 
half occupy the O sites together with the A cations. In 
general, spinels can be characterised by the degree of in¬ 
version (i.e. the concentration of B cations occupying T 
sites). 

The bulk structure (unstrained case) of the inverse 
spinel CoFe 2 0 4 , with a disordered Co 2+ /Fe 3+ occu¬ 


pancy on the octahedrally coordinated sites, has cu¬ 
bic symmetry (space group Fd3m). This enforces the 
quadratic magnetization terms in the magnetocrystalline 
anisotropy to vanish, leaving as lowest order terms the 
quartic ones, which are typically of smaller size, although 
they can be appreciably large 11 ^ However, strain lifts 
the cubic symmetry, allowing quadratic magneto crys¬ 
talline anisotropy, which can either favor in-plane or out- 
of-plane magnetization. While in CoFe 2 0 4 compressive 
strain favors in-plane magnetization, i.e. magneto c ry s- 
talline and shape anisotropies cooperate 7 1Q I 12 I 17 I in 
CoFe 2 0 4 films grown under tension on MgO, the strain- 
induced out-of-plane bulk anisotropy is large enough to 
overcome the shape anisotropy, promoting out-of-plane 
magnetization^ 9 13 . The magnetocrystalline anisotropy 
in CoFe 2 0 4 is expected to be th e re sult of strong spin 
orbit coupling at the Co 2+ cation 18 19 and, thus, a signif¬ 
icant effect is expected for different site occupation of the 
Co 2 + lonP® as observed in CoFe 2 0 4 under different 
annealing conditions 10 . 

Another Co 2+ compound with cubic spinel structure, 
CoCr 2 0 4 , raised considerable interest a few years ago 
as it was discovered to be a multiferroic with strong 
magnetoelectric coupling when unde rgoing a transition 
to a ferrimagnetic conical spiral state^H 25 . As temper¬ 
ature is lowered, CoCr 2 0 4 first shows a transition to a 
collinear ferrimagnetic state at Tjvi~ 95K, which exists 
down to TN 2 r ^ 1 26K. Below this temperature the ferri¬ 
magnetic conica l spir al state appears and induces an elec¬ 
tric polarization 22 25 . At T^ 14K an abrupt change in 
electric polarization occurs 22 The fact that it is diffi¬ 
cult to grow good quality single crystals l eaves some open 
questions on the detailed magnetic stat d 25 | 26 J . To date, 
the epitaxial growth of CoCr 2 0 4 in thin film form has 
only been reported by a few groups 27 30 . No study so far 
has reported strain engineering of magnetic properties. 

CoCr 2 0 4 has a strong tendency to order in a normal 
spinel structure 3 ^, i.e. with Co 2+ occupying T sites and 

Cr 3 + 

occupying the O sites. This is different from the 
behavior of CoFe 2 0 4 and NiFe 2 0 4 , which have an in- 
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verse spinel structure, with the Co 2+ and Ni 2+ occupy¬ 
ing O sites, while Fe 3+ cations are distributed in equal 
proportion among O and T sites (although real samples 
can show partial occupation of Co 2+ in T sites due to 
growth kinetic effects®). Because of the different Co 2+ 
distribution in CoCr 204 and CoFe 204 , the effect of strain 
on magnetic anisotropy is expected to be different in the 
two compounds. However, it is unknown how these differ¬ 
ences will manifest themselves in actual films, since stud¬ 
ies on the magnetic anisotropy of strained CoC^CU, as 
well as a comprehensive theoretical treatment for spinels, 
are lacking. 

In this work, we investigate experimentally and the¬ 
oretically the strain dependence of magnetocrystalline 
anisotropy in CoCr 204 , with the purpose of generalizing 
the behavior of spinel systems in which Co 2+ is expected 
to be the main source of magnetocrystalline anisotropy. 
After a brief discussion of experimental and computa¬ 
tional details, we present our results of magnetic mea¬ 
surements performed for thin films of CoC^CU under 
both compressive and tensile strain. In particular, we 
find that in CoC^CU strain affects magnetocrystalline 
anisotropy in the opposite way compared to CoFe 204 i 
compressive and tensile strain favor, respectively, out-of- 
plane and in-plane magnetization. We then discuss re¬ 
sults of ab initio calculations of the strain dependence of 
magnetocrystalline anisotropy. Finally, we provide sim¬ 
ple arguments, based on second order perturbation the¬ 
ory, which explain the opposite tendencies of the strain 
dependence of magnetocrystalline anisotropy in CoCr 204 
and CoFe 204 . 


I. EXPERIMENTAL AND COMPUTATIONAL 
DETAILS 

CoCr 2 04 thin films were grown using pulsed laser de¬ 
position (PLD). The RHEED (Reflection High-Energy 
Electron Diffraction)-assisted PLD system is equipped 
with a Lambda Physik COMPex Pro 205 KrF laser with 
a wavelength of 248 nm. The laser was run with a fre¬ 
quency of 0.5 Hz for 2 hours, using a laser fluence of 3 
J/cm 2 and a target-substrate distance of 50 mm, which 
resulted in CoC^CU films with thicknesses of about 
40 nm. The substrate was kept at 400-450 °C during 
growth. A ceramic CoCr 204 pellet, sintered by solid 
state reaction^, was used as target. The growth took 
place in a 0.01 mbar oxygen plasma atmosphere, created 
by an Oxford Scientific MPS-ECR plasma source in order 
to improve the oxidation of the films 33 . After deposition, 
the films were cooled down in 0.5 bar O 2 atmosphere with 
a cooling rate of 5 °C per minute. The structure of the 
thin films was determined by X-ray diffraction (XRD), 
using an X’pert lab diffractometer. Magnetization was 
measured using a QD MPMS 7 XL SQUID magnetome¬ 
ter. 

Thin films were deposited both on MgO and MgA^CU 
substrates. The substrates were post-annealed at 650 °C 
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FIG. 1 . (color online) (a) Crystallographic structure of 
CoCr204- The primitive unit cell of the fee lattice is shown, 
containing 2 formula units. Blue, green, and red spheres rep¬ 
resent, respectively, Co 2+ , Cr 3+ and O 2- ions. Grey arrows 
indicate the ferrimagnetic ordering considered in the calcu¬ 
lations. Red arrows in panels (b) and (c) show the oxygen 
displacements for the two modes which appear on relaxing 
the structure under strain. 


overnight 34 in lxlO -5 bar oxygen atmosphere. MgAUCU 
has the spinel crystal structure, like CoC^CU, but with 
a smaller lattice parameter, corresponding to a lattice 
mismatch!^ of 3.1 % that induces a large in-plane com¬ 
pressive strain in the CoC^CU films. MgO has the rock- 
salt structure with a lattice mismatch of 1.1 %, resulting 
in CoCr 204 films under tensile strain. 

Ab initio calculations were performed using the 
LSDA+U approach, as implemented in the Vienna Ab 
initio Simulation Package (VASP), with Projector Aug¬ 
mented Wave pseudopotentials. We use a 7 x 7 x 7 
Monkhorst k-point grid and consider the primitive unit 
cell of the underlying face-centered cubic (fee) lattice, 
containing two formula units. The values of the effective 
on-site Coulomb interactions and Hund’s couplings for 
Co 2+ and Cr 3+ were set to Uc 0 = 4 eV, Jq g = 1 eV 
and Ucr = 3 eV, J§ r = 1 eV, respectively. As shown in 
Ref. [313 a realistic description of the magnetic properties 
of CoCr 2 04 is obtained using these values. Moreover, 
for both transition metal ions, the 3 s and 3 p states were 
treated as valence electrons. The lattice constants and 
atomic positions, both for the unstrained cubic unit cell 
and for various values of strain, were obtained by calcu¬ 
lating the energy of relaxed structures at different values 
of the out-of-plane lattice constant and extracting the 
value corresponding to minimal energy. Relaxations were 
performed considering a collinear ferrimagnetic state in 
which the spins of the Co 2+ sublattices are antiparallel 
to those of the Cr 3+ sublattices (see Fig. Eh 

Calculations to obtain the size of the magnetocrys- 
talline anisotropy were performed using non-collinear 
spin-polarized calculations, in which the direction of the 
spins of each magnetic sublattice was constrained and ro¬ 
tated, keeping the ferrimagnetic order described above. 
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FIG. 2. (color online) (a) Specular XRD pattern around 
the (004) CoCr 204 reflections. The inset shows the rock¬ 
ing curves around the (004) reflections for both films. XRD 
Reciprocal space maps (RSMs) around the (115) reflections 
of the CoCr 204 films grown on MgADCU (b) and MgO (c). 
The axes in b) and c) are labelled in reciprocal lattice units 
(r.l.u.) of the substrate. The vertical dashed line in c) marks 
the position of the substrate peak (not shown) 


II. RESULTS AND DISCUSSION 
A. Experimental results 

Despite the lattice mismatch between both substrates 
and CoCr 204 , high quality epitaxial thin films have 
been grown under both tensile strain, on [001]-oriented 
MgO substrates, and under compressive strain, on [001]- 
oriented MgAl 204 substrates. The possibility of apply¬ 
ing high amounts of strain to spinel materials was also 
observed in CoFe 2 04 by Gatel et a/. 12 . The XRD charac¬ 
terization of the films shows a large number of thickness 
fringes around the film peaks, revealing a superb flatness 
of the film interfaces (Fig. IT)- The excellent crystalline 
quality of the films is shown by the rocking curves of 
the film peak, displaying full width at half maximum 
(FWHM) below 0.03° (inset of Fig. [2^). The films on 
MgAl 2 04 and MgO have out-of-plane lattice constants of 
8.58 A and 8.17 A, respectively (bulk value is c = 8.33 A). 

The reciprocal space map around the (115)-reflection 
of CoCr 2 04 in Fig. [ 2)3 shows that the in-plane lattice con¬ 
stant of the films are identical to that of the MgA^CU 
substrate, as substrate and film peaks are observed at 
the same value of qHKO • Thus, the films are coherently 
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FIG. 3. (color online) Magnetization of strained CoCr 204 on 
MgO and MgAl 204 , measured both in-plane and out-of-plane, 
in a field of 0.01 T after cooling in 1 T. 


strained, despite the huge lattice mismatch between the 
MgAl 2 04 and bulk CoCr 2 04 of 3.1%. Despite the rel¬ 
atively low intensity of this reflection, thickness fringes 
are also visible around the (115)-film peak. 

According to our specular XRD measurements, 
CoCr 204 films grown under tensile strain on MgO dis¬ 
play a similar crystal quality to that of the films grown 
on MgAl 204 . The (115)-off-specular peak of CoCr 204 in 
Fig# shows that the film on MgO is also fully strained. 
However, the rocking curve of this peak [scan in the qnHO 
direction], has a FWHM of 1.2°. This is the result of the 
difference in structure between the rock-salt MgO sub¬ 
strate and the spinel CoCr 2 04 film. It is known that 
the growth of spinel films on rock-salt structure materi¬ 
als can result in the formation of anti-phase boundaries, 
creating a less homogeneous lattice in the in-plane direc¬ 
tions. Antiphase boundaries have, indeed, been observed 
in our CoCr 204 thin films using TEM. Their formation 
and effects will be reported elsewhere. 

Magnetization of the thin films was measured in both 
in-plane and out-of-plane geometries. The magnetiza¬ 
tion of CoCr 204 is highly anisotropic as shown in Fig. [3] 
Films on MgA^CU, under compressive strain, in a field 
of 0.01 T applied perpendicular to the film, show the 
ferrimagnetic transition at 81 K (the bulk value be¬ 
ing T/v 1 ^ 95 K). Applying a magnetic field along the 
[100] in-plane direction induces only a paramagnetic tail, 
which originates from the substrate contribution, show¬ 
ing that the film has an easy out-of-plane axis. The films 
on MgO, under tensile strain, show a ferrimagnetic tran¬ 
sition at ^80 K under a 0.01 T field applied in the [100] 
in-plane direction, while no ferrimagnetic transition is ob¬ 
served when the field is applied in the out-of-plane direc¬ 
tion. The magnetic anisotropy in CoCr 204 films, thus, 
changes sign depending on the sign of the strain: the 
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FIG. 4. (color online) Magnetic hysteresis loops of CoC^CP 
on MgAl 2 04 , under compressive strain (a), and on MgO, un¬ 
der tensile strain (b), at 20 K, measured in the (100) in-plane 
and (001) out-of-plane directions. 

magnetic easy axis is out-of-plane in films under com¬ 
pressive strain, while it is in-plane in films under ten¬ 
sile strain. As discussed above, the strain-tuning of the 
magnetic anisotropy has also been observed in films of 
CoFe 204 but in that case the effect is opposite 7 . 

The magnitude of the magnetization is also dependent 
on the sign of the strain. The low-temperature magne¬ 
tization of the CoCr 2 04 films under compressive strain 
is 0.34 /Ib per formula unit (f.u.), while it is only 0.03 
11b/ f.u. under tensile strain. These values are, respec¬ 
tively, larger and smaller than the 0.08 /i#/f.u. found 
in single crystald^. Cation inversion is well known to 
increase or decrease magn etizatio n in spinel structures 
as for instance in NiFe 2 O j 20 * 37 * 3 ^. However, NiFe 204 
and CoFe 204 have a collinear ferrimagnetic structure, 
whereas in bulk CoCr 204 , the low temperature mag¬ 
netic structure is a ferrimagnetic conical spiral in which 
the low temperature magnetization is not determined by 
cation inversion (Co 2+ and Cr 3+ both have three un¬ 
paired spins) but by the conicities of the magnetic spirals 
of the various sublattices. 

A ferrimagnetic conical spiral can also be present in 


thin film samples at low temperatures. This, together 
with a possible dependence of the spiral conicities on the 
applied strain, might explain the change in the low tem¬ 
perature magnetization. It is also worth to note that 
the gyromagnetic tensor is likely to depend weakly on 
strain, which might contribute to the discrepancy be¬ 
tween the low temperature magnetization of the two 
samples. Moreover, the magnetization of spinel struc¬ 
tures on MgO can be smaller than in bulk as a result of 
anti-phase boundaries 9 39 40 . As mentioned above, anti¬ 
phase boundaries were observed in our CoCr 204 films on 
MgO, and could thus explain their lower magnetization 
values. However, since no APB have been observed in our 
films on MgAl 2 04 , it is likely that the increased value of 
magnetization with respect to bulk CoCr 204 in this case 
has other origins, such as changes in conicities. 

The strain-dependent magnetic anisotropy is con¬ 
firmed by the magnetic hysteresis measurements shown 
in Fig. [4j Films on MgAl 204 (Fig. @0 measured with 
the field along the out-of-plane direction, show a square 
loop with coercive field of 2.2 T when the diamagnetic 
slope of the substrate contribution is subtracted, while in 
the in-plane direction no loop opens. This confirms that 
the [001] out-of-plane and [100] in-plane axes are, indeed, 
the easy and hard axes, respectively. The hysteresis loop 
on MgO has an opening when the field is applied in the 
[100] in-plane direction while no loop opening is observed 
in the out-of-plane direction (Fig. [ 4 ) 3 ). On MgO, the 
loop in the in-plane direction is not square, which could 
be a sign of the field not being applied along the easy 
axis. However, it is well known that films of CoFe^O^ 
with anti-phase boundaries have such loop shapes 17 39 . 
In our CoCr 2 04 films on MgAl 204 , we observe an ex¬ 
tremely large coercive field (unusual for spinels), which 
most likely originates in the symmetry lowering induced 
by the strain. Extremely high coercive fields in spinel 
thin films have also been observed by Yanagihara et alP 1 . 


B. Ab initio calculation of magnetic anisotropy 

To investigate the microscopic mechanisms leading to 
the strain dependence of the magnetic anisotropy shown 
in Fig. [3j we perform ab initio calculations. We first 
calculate the crystal structure under strain by relaxing 
the out-of-plane lattice constant together with the atomic 
positions to their lowest energy value for a set of in-plane 
lattice constants. We then use the relaxed structures to 
calculate the magnetocrystalline anisotropy. 

In addition to the change of the lattice constants, the 
structural relaxation of the strained unit cell induces two 
structural distortions involving oxygen displacements. 
These correspond to displacement modes of symmetry 
T+ and T+. Mode corresponds to a uniform com¬ 
pression (enlargement) of the oxygen tetrahedra coordi¬ 
nating Co 2+ when compressive (tensile) strain is applied 
(see Fig. [I]d). Mode r+, on the other hand, changes the 
local symmetry from that of a perfect tetrahedron to that 
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FIG. 5. (color online) The left and right panels show, respec¬ 
tively, the calculated angle ( 0 ) dependence of the energy (with 
the energy at 0 — 90° subtracted) for —3.12% (compressive) 
and 1.08% (tensile) strain. 6 — 0 corresponds to magnetic mo¬ 
ment along the x axis (in-plane), while 0 m 90° corresponds to 
magnetic moment along z (out-of-plane). The two values of 
the strain are approximately equal to the experimental ones. 

of a tetragonally distorted one, as shown in Fig. 0 )- We 
note that the displacements of O 2- in mode give rise 
to shifts in the crystal field splitting similar to those cor¬ 
responding to the change of the lattice constants when 
strain is applied. For example, mode when com¬ 
pressive (tensile) strain is applied, shrinks (enlarges) the 
tetrahedra surrounding Co 2+ in the xy plane and elon¬ 
gates (shrinks) them along the z axis. 

To estimate magnetocrystalline anisotropies we con¬ 
sider the anisotropic part of the magnetic Hamiltonian 

h a = Y j ^ ) (Si v ) 2 + 

£ A ( 4)((^) 4 + (^ !/ ) 4 + (^) 4 ), ( l ) 

where, i and v = 1,6 label, respectively, the unit 
cell and the magnetic sublattices in the unit cell, A ^ 
and A! 4 ) are, respectively, the strength of quadratic and 
quartic anisotropy of spins at sublattice v and S a , with 
a = x,y,z, are the spin components along the pseudocu- 
bic crystallographic axes. Here we denote with z the 
direction perpendicular to the film. On the one hand, 
for the bulk case, the first term on the right-hand side 
of Eq. ([l]) has to vanish due to the cubic symmetry. On 
the other hand, while the fourth-order term is always 
allowed, relatively small tetragonal distortions are ex¬ 
pected to make the second-order term dominating. 

We map ab initio energies to the Hamiltonian in Eq. 0 
using the following procedure (see e.g. Ref. 02|. We per¬ 
form a set of ab initio calculations where the direction of 
all the spins is rotated uniformly (i.e. keeping the mag¬ 
netic ordering of Fig. [l]) from the in-plane direction {x) 
to the out-of-plane direction (z) by an angle 6. For these 
magnetic configurations, Eq. (HI) yields a 0 dependence 
of the ab initio energy per unit cell given by 

e{6) = a cos (0) 2 + /? (cos (#) 4 + sin(#) 4 ) + 7 . ( 2 ) 

Here, 7 is a constant energy shift which includes all con¬ 
tributions unrelated to magnetocrystalline anisotropy, 
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FIG. 6. (color online) Calculated strain dependence of the 
constant o, as described in the text. The inset shows the 
strain dependence of the constant f3 related to the quartic 
anisotropy. 


while a and (3 include, respectively, the quadratic and 
quartic contributions to spin anisotropy of all six mag¬ 
netic sublattices in the unit cell. 

The two panels in Fig. [5] show the numerical results 
(dots) and the fits (curves) obtained using Eq.pl) for e(0) 
for strains close to the experimental values of —3.1 % 
(compressive) and 1.1% (tensile). We note that e(0) fits 
well the ab initio energies, indicating that the model in 
Eq. 0 describes well the magnetocrystalline anisotropies 
in this compound. Moreover, we find that, in agreement 
with experimental results, a tensile value of the strain 
favors a direction of the spins in the film plane (energy 
minimum at 0 = 0°), while compressive strain favors a 
direction perpendicular to the film plane (energy mini¬ 
mum at 0 = 90°). 

For the case of a collinear ferrimagnetic state, the mag¬ 
netocrystalline anisotropy part of the microscopic mag¬ 
netic Hamiltonian, Eq. 0, can be easily related to the 
direction of magnetization, i.e. magnetization is favoured 
along the easy axis anisotropy direction. However, for a 
ferrimagnetic conical spiral state, such straightforward 
relationship does not always hold. In this state, the 
magnetic ordering of the various sublattices consists of a 
collinear and a spiral component and the collinear com¬ 
ponent, which gives rise to the net magnetization, is per¬ 
pendicular to the plane containing the spiral one. There¬ 
fore, whether the collinear or the spiral part of the mag¬ 
netic ordering is along the easy axis direction might be 
related to the conicity of the spin ordering, i.e. the rel¬ 
ative size of the collinear component with respect to the 
spiral one. Nonetheless, we note that the transition tem¬ 
peratures in Fig. [3] are comparable with the transition 
temperature of the collinear ferrimagnetic state in the 
bulk (T/vi ~ 95K). It is thus reasonable to assume that, 
at least at temperatures below but close to the transi¬ 
tion, the magnetic ordering is collinear also for the films 
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and, thus, the magnetization is parallel to the micro¬ 
scopic easy direction given by Eq. 0 - 

To generalize the dependence of the magnetocrystalline 
anisotropy on strain, we extract the values of the con¬ 
stants a and /? as the strain is tuned. These dependences 
are shown in Fig. [6j We note that the strength of the 
quadratic part (a) of the magnetocrystalline anisotropy 
decreases linearly (to first approximation) as the strain 
is tuned from compressive to tensile and changes sign for 
the bulk structure. In contrast, the strength of the quar- 
tic term (/3) shows a much weaker dependence on strain. 

It is reasonable to assume that the strongest contribu¬ 
tion to the magnetocrystalline anisotropy originates from 
the Co sublattice due to the d 7 electronic configuration of 
Co 2+ and the larger value of spin-orbit coupling (which 
increases with the atomic number) compared to Cr 3+ . To 
verify this assumption, we estimate the separate contri¬ 
butions of the Co 2+ and Cr 3+ sublattices by performing 
the following calculations. We constrain the spins of the 
Cr sublattices to be ferromagnetically aligned along the 
y direction and the spins of the Co sublattice to be fer¬ 
romagnetically aligned along a direction in the xz plane 
forming an angle Oqo with the x axes. As Oqo is changed, 
the exchange energy within the Co (Cr) sublattices can¬ 
not vary, as all the spins within these sublattices are kept 
parallel. Moreover, the exchange energy between the Co 
and Cr sublattices is always zero, as their spins are kept 
perpendicular to each other. Therefore, the main contri¬ 
bution to the spin part of the Hamiltonian should origi¬ 
nate from Eq. 0 , which, for such a spin configuration, 
has the energy dependence 


e(0Co) = «co cos(#co) 2 + /3co (cos(0 C o) 4 + sin(# Co ) 4 ) • 

(3) 

Here, ac 0 and /3c 0 contain the contribution to the 
anisotropy of the Co 2+ sublattice and a constant energy 
shift is omitted. Similar calculations can be performed 
to obtain the anisotropy contributions of only the Cr 3+ 
sublattice, ac T and /3c r - 

Figure[7](a) shows the comparison of the contributions 
of the Co and Cr sublattices to the magnetocrystalline 
anisotropy for the case of 3% tensile strain. In agree¬ 
ment with our initial assumption, the magnetocrystalline 
anisotropy is dominated by the contribution of Co ions. 
Interestingly, we note that the sum of the contributions 
of the two sublattices does not give exactly the value ob¬ 
tained for the total magnetocrystalline anisotropy. We 
further investigated the origin of such difference and ex¬ 
cluded possible sizeable contributions from terms such as 
Dzyaloshinskii-Moriya interaction or anisotropies of the 
exchange interactions. A possible source of this discrep¬ 
ancy might arise from a non-local contribution to the 
magnetocrystalline anisotropy, e.g. a contribution re¬ 
lated to hopping between Co 2+ and Cr 3+ which might 
depend on the relative orientation of the spins of the two 
sublattices. 



FIG. 7. (color online) (a) Calculated energy dependence on 
the angle 0 when the contributions of the Cr 3+ (green squares) 
and Co 2+ (light blue diamonds) are separated as described 
in the text for 3% tensile strain. The sum of the two con¬ 
tributions and the energy dependence obtained by a global 
rotation of all the spins are shown by (violet) triangles and 
(blue) dots, respectively. Curves represent the fits obtained 
as described in the text, (b) Crystal field splitting of a Co 2+ 
ion in tetrahedral coordination under compressive strain. 


C. Second-order perturbation theory for the 
magnetocrystalline anisotropy 


After establishing that the dominant contribution to 
magnetocrystalline anisotropy originates from the Co 2+ 
sublattices, we note that the strain dependence of the 
anisotropy can be understood using a simple picture of 
isolated Co 2+ ions coordinated by distorted oxygen tetra- 
hedra. In this case the Co 2+ contribution to the first term 
on the right-hand side of Eq. 0 can be obtained by sec¬ 
ond order pert urbation theory in the spin-orbit coupling 
constant ^ 43 | 44 l This gives the expression, 

H a = & { cl{S z ) 2 = -\ 2 (\ zz - A XX )(S Z ) 2 , (4) 

where we treat the Co spin as a classical vector and 
A n = - • Here, Li is the i-th component of the 

orbital momentum operator, | k) and |0) denote, respec¬ 
tively, the k- th excited state and the ground state. To ob¬ 
tain Ha , we use the tetragonal symmetry, i.e. A xx = A yy . 
We first note that, for the case of a (“normal”) cubic 
spinel the symmetry of the ions sitting in the T sites 
is Trf. Considering the joint effect of the change of the 
lattice constants and the mode Tg", one obtains (for the 
case of compressive strain) the splitting of the d levels as 
sketched in Fig. 0(b). Here, we consider the d 7 electronic 
configuration and we assume that A > A t 2 > A e > 0, as 
the lobes of £2 states are pointing more directly towards 
the ligands. We also note that, in principle, crystal field 
splittings can be estimated from the ab initio calculations 
using Wannier functions. However, it was shown (see e.g. 
Ref. 05) that this method requires a careful analysis and 
the so-obtained splittings do not necessarily correspond 
to a simple point charge model of isolated ions within a 
crystal field. By taking into account only the electronic 
excitations with lowest energy, i.e. only one-electron ex- 
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citations not violating Hund’s rule, we obtain 

a ( 2 ) _ A 2 [—4A t2 (A + A t2 ) + 3A e (A + A e )] 

Co ' T _ (A + A t2 ) (A + A e ) (A + A e + A t2 ) ’ 

where the subscript T indicates that the equation holds 
for Co 2+ in tetrahedral coordination. This expression, for 
the case of compressive strain, leads to an out-of-plane 
easy axis anisotropy (AJ^ is negative as A t 2 > A e , which 
favors spin orientation along z), in qualitative agreement 
both with experiments and ab initio calculations. More¬ 
over, we note that Eq. © holds also for the case of tensile 
strain, where A t 2 and A e are expected to be negative 
with A t 2 < A e < 0. This gives rise to an easy-plane 
value of the quadratic anisotropy (positive A^, which 
favors spin orientation perpendicular to z). 

Furthermore, we speculate that a similar argument 
can be used to explain the sign of the magnetostriction 
in the inverse spinel CoFe 204 . Assuming that also in 
CoFe 204 the largest contribution to the magnetocrys- 
talline anisotropy originates from the Co 2+ ions, we 
use similar considerations to those used for the tetra- 
hedrally coordinated Co 2+ to obtain an expression for 
the anisotropy for the case of octahedral coordination. It 
is important to note that in cubic spinel (both normal 
and inverse) the octahedral sites do not possess symme¬ 
try O h but are trigonally distorted. Such distortion was 
previously shown to play an important role in the surpris¬ 
ingly large quartic magnetic anisotropy found in cubic 
Co a ,Fe 3 _ ;E 04 18 19 . Nonetheless, as we are interested in 
the emergence of second order anisotropy with a strain- 
induced tetragonal distortion, we consider the case of a 
perfect O h symmetry for the unstrained case. Moreover, 
we restrict ourselves to the case of compressive strain 
where, similarly to the case of CoC^CU, non-degenerate 
perturbation theory can be applied^. For compressive 
strain, we obtain 

A (2) _ A 2 [A t 2 g (4A eg + 5A t2 g) + A (A + A eg + 6A t 2 g)\ 

C ° ° A^2g (A + A t 2 g) (A + A eg + A t 2 g) 

( 6 ) 

where the subscript O indicates octahedral coordination, 

A t 2 g &xy £xz/yzi A eg S x 2_ g 2 £ z 2 , A S z 2 &xy 

and £i is the on-site energy of the i-th d orbital. For 
compressive strain of the octahedron, it is reasonable to 
assume that A eg > 0, A t 2 g > 0 and A > 0, which implies 


that the anisotropy, Eq. ([6]), favors in-plane orientation 
of the spins (A^_ Q > 0). We stress that in contrast to 
Eq. the expression in Eq. © does not hold for tensile 
strain as the treatment of this case would require the use 
of degenerate perturbation theory. 

III. SUMMARY 

We report the growth of fully strained (001) CoCr 2 04 
films on MgA^C^ and MgO substrates. Magnetometry 
experiments reveal that compressively strained CoCr 204 
thin films exhibit a perpendicular magnetization, while 
films under tensile strain exhibit an in-plane magneti¬ 
zation. The mapping of LSDA+U calculations to the 
anisotropic part of an effective spin model, Eq.Q, shows 
that (001) compressive strain favours a spin easy axis 
along the [001] direction while tensile strain favours spin 
orientation in the film plane. The quadratic anisotropy 
term changes smoothly as strain is tuned and vanishes in 
the cubic case. Moreover, these calculations confirm that 
the main contribution to magnetocrystalline anisotropy 
is given by the Co 2+ sublattice. The direction of magne¬ 
tization is parallel to the easy direction of the microscopic 
magnetic model for the collinear ferrimagnetic state and 
for ferrimagnetic conical states with large enough conic- 
ity. We note that, for both the samples, the ferrimag¬ 
netic state appears at temperature close to T^i below 
which bulk CoCr 2 04 is collinear for a large temperature 
regime (26K < T < 95 K). Therefore, it is reasonable 
to assume that the magnetization is parallel to the easy 
direction of model Eq. (fTh . Under this assumption the 
results obtained with LSDA+U are in agreement with 
experimental observations. 

Furthermore, to explain the behavior of the quadratic 
part of the magnetocrystalline anisotropy with strain in 
CoCr 2 04 , we give a simple argument based on a pertur¬ 
bative treatment of spin-orbit coupling when the crys¬ 
tal field splitting is modified by strain. Finally, we use 
these arguments to explain why the strain dependence of 
anisotropy is opposite in the inverse spinel CoFe 204 . The 
same sign of the strain-induced anisotropy as in CoFe 204 
has been reported in CoO filnd^U (also containing Co 2+ 
in octahedral coordination), indicating that the simple 
argument given in our discussion might not be only spe¬ 
cific for spinel structures but can, in principle, be ex¬ 
tended to more general cases. 
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